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Q ■ Abstract 

I Stochastic variational inequalities provide a unified treatment for stocliastic diS'eren- 

0^ ' tial equations living in a closed domain with normal reflection and (or) singular repellent 

drift. When the domain is a polyhedron, we prove that the reflected-repelled Brownian 
' motion does not hit the non-smooth part of the boundary. A sufficient condition for non- 

p ^ I hitting a face of the polyhedron is derived from the one-dimensional case. A complete 

answer to the question of attainability of the walls of the Weyl chamber may be given for 
a radial Dunkl process. 



1 Introduction 



^ ' There have been many works about stochastic differential equations with reflection on the 

. boundary of a domain. In some of them the domain is a convex polyhedron (|17j. [29], |30j . 

|12j]. A typical question in this setting is the following: does the continuous process 
hit the non-smooth part of the boundary? The answer depends on the drift and diffusion 
coefficients of the process and on the direction of reflection (normal or oblique). In particular, 
R.Williams ^30j has proven that the Brownian motion with a skew symmetry condition on 
the direction of reflection does not touch the intersections of the faces of the polyhedron. 

On the other hand there also exists an extensive literature about non-colliding Brownian 
particles ( [15] , [3] , [18] , [16] , [23] ) . Most of these works originate in the study of the eigenvalues 
^ ' of Gaussian matrix processes. These eigenvalues are solutions to systems of stochastic differ- 

■ ential equations with a singular drift that prevents the particles from colliding. Extensions 

of these systems are Dunkl processes |25j that have recently been developed in connection 
with harmonic analysis on symmetric spaces. The radial part of a Dunkl process may be 
considered as a Brownian motion perturbed by a singular drift which forces the process to 
live in a cone generated by the intersection of a finite set of half-spaces ( [9] , |10j ) . Depending 
on the values of some parameter, the process may touch the walls of the cone or not. 

Actually it is possible to unify both theories of (normal) reflection and strong repulsion 
within a common framework. This is the role of stochastic variational inequalities, also 
called multivalued stochastic differential equations (MSDE) that were mainly developed by 
E.Cepa ([4], [5]). These equations are associated to a convex function in a domain of R'^. 
Depending on the boundary behavior of this function the diffusion will (normally) reflect on 
the boundary, hit the boundary without local time, or live in the open domain. We shall 
here follow this way and concentrate on a Brownian motion living in a convex polyhedral 
domain, bounded or unbounded. To each face of the polyhedron is associated a repelling 
force with normal reflection when the repulsion is not strong enough. In this setting we shall 
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ask whether the process may hit the various faces. Our first task will be to rule out the 
possibility of hitting the intersection of two faces. Once this is achieved, the problem is now 
basically one-dimensional and we may use the ordinary scale function of real diffusions. 

In several previous works ([20j, [8j), this issue has been studied in the particular case of the 
hyperplanes Hij := {x = (xi, . . . , Xd) G IR"^ : Xi = Xj}, i j and presented as the problem of 
collisions between Brownian particles. There is a simple collision if two coordinates coincide 
and a multiple collision if at least three coordinates coincide at the same time. Because 
the d-dimensional Brownian motion does not hit the intersection of two hyperplanes, one 
can guess that an additional drift does not change anything. However a rigorous proof is 
necessary because the singularity of the drift makes useless the usual Girsanov change of 
probability measure. The counterexample of Bass and Pardoux [T] also showed that uniform 
nondegeneracy of the diffusion term does not preclude multiple collisions. 

As in [8] where the particular case of electrostatic repulsion was considered, our proof only 
uses basic tools from stochastic calculus, mainly McKean's martingale method [22] which was 
already used in [2\ to prove non-collision for the eigenvalues of Wishart processes. Another 
way could be to use the theory of Dirichlet forms as done in ^20j where a general condition 
of non-collision has been obtained. 

The paper is organized as follows. In Section 2 we introduce basic definitions and nota- 
tions. The main features about stochastic variational inequalities are also recalled. Section 
3 is devoted to non attainability of the edges of the polyhedron. In Section 4 we give a suffi- 
cient condition of non attainability of a single face. Section 5 presents some applications to 
Brownian particles with nearest neighbor interaction, Wishart processes and Dunkl processes. 



2 Multivalued stochastic differential equation in a polyhedral 
domain 

Let {n,J^,{J-'t,t > 0),P) be a filtered probability space endowed with the usual conditions 
and B = {Bt) be a (.7-i)-adapted d-dimensional Brownian motion starting from the origin. 
Let 

$ : M'^ ^ (-oo,+oo] (1) 
be a lower semi-continuous convex function such that 

dom{(^) := {x : ^{x) < oo} (2) 

has nonempty interior. Let 

D := Int{dom{<^)) . (3) 

For simplicity of notation, we will assume that $ is on D. If x G dD, we say that the 
unit vector n{x) is a unit inward normal to at x if 

n{x).{x - z) <0 (4) 

for any z € D. Based on the results in [4j, the following theorem has been proved in [6] (see 
also Theorem 2.2 in [7]). 

Theorem 1 For any J^Q-measurable random variable Xq with values in D, there exist a 
unique continuous [Tt)- (adapted process X = {Xt,Q < t < oo} with values in D and a unique 
continuous {Tt)-o-dapted non- decreasing process L = {Lt,0 <t< oo} such that 



Xt = Xq + Bt - /o V$(X,) ds + /o UsdLs < t < oo 

= Jo '^{x,edD}dLs < t < oo 



(5) 
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where Ug is dLg-a.e. a unit inward normal to D at Xg- For any <T < oo, 

/ ^{XsedD}ds = (6) 

^0 



and 



[ \V'^{Xs)\ds < oo . (7) 
Jo 



From now on we concentrate on a particular polyhedral setting. Let I := {1, . . . ,m} 
where m > 1. We consider a convex function $ of the following form 



<^{x) := J2Mx-ni - ai) (8) 
iei 



where for any i €l, 



(pi is a convex l.s.c. function , (pi = +oo on (— oo,0), (pi is on (0,+oo) 

Hi is a unit vector (9) 

Ui is a real number . 

We may assume all rii are different. Then, 

V^{x) = J2iei'^i4''i{x.ni - ai) 

D = {x G M'^ : x.Ui > ai yi£ 1} (10) 

D = {x G M°' : x.Ui > ai Vi G 1} . 

As D is not empty, there exists a ball with center y £ D and radius b > included in D. Let 
Xf be the solution given by Theorem 1. For i G I let 

Ul := Xt.m - ai . (11) 

We will need a strengthening of inequality (7) ([7J,Th.2.2). 

Lemma 2 For any i £ I, for any < t < oo, 



t 



{Ul)\ds < oo . (12) 







Proof. This is clear if (p'-{0+) > — oo. Let 

J:={jGl:(/>;.(0+) = -oo} (13) 
and let < e < 6 be such that (p'j{u) < for any j G J and u G (0, e). For K C J let 

Ak ■■= {x eD : x.Uj < aj +£ Vj G K, x.Uj > Oj + e Vj G J \ K} . (14) 
Then for t > 

r IaA^s)] n, (P'^mids <Y. f UA^sM'^mids < oo . (15) 



3 



Using (7) we get 

/ E 'f^'jiUDlds < oo (16) 

•''0 



and therefore 



< JiiAAXs)\y-x,\^^^^\<f>'.{ui)\ds (17) 

< oo 

from the continuity of X on [0,t\. Then for any j G J 

j^w,{ui)\ds = /o*v<,}i</';(f/^)id. + /*v>,jK(f/i)id. 

= EjeKcjIo^AAXs)\<^^{m)\ds + J^l^^,^^^^\c^^{Ui)\ds (18) 
< oo . □ 

For any J C I, J 7^ 0, we set 



Hj 



= {x eW^ : x.rij = aj Vj G J} 

= {x G M!^ : x.rij = aj Vj G J, x.nj > aj Vj J} 

= inf{t >0 : Xt e Hj} 

= mi{t >Q : Xt eKj} . 



(19) 



Lemma 3 Let J C I and V := span{nj,j G J}. // n(x) is a unit inward normal to D at 
X G Kj, then n{x) G V. 



Proof. Let v J-V. For e > small enough, 

zi = X + ev Z2 = X — ev 



satisfy 



Then 



and therefore 



zi.Uj = aj Vj G J zi.Ui > a-i Vi J 
Z2-nj = aj Vj G J Z2-ni > ai Mi ^ i 

n{x).{x — zi) < n{x).{x — -22) < 

n{x).v = . 



□ 



3 Nonattainability of the edges 

This section is devoted to the proof of the following theorem. 
Theorem 4 For any J C I with card{3) > 2, 

P((7J = 00) = 1 . 
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Proof, a/ We first consider the initial condition Xq. Prom (6) we deduce that for any u > 

there exists < v < u such that X,, G D a.s. Using the continuity of paths and the Markov 
property we may and do assume that Xq G D in order to prove that aj = oo a.s. 
b/ We will also assume that 

max(/>^(0+) < . (20) 

If not we introduce for any < T < oo the equivalent probability measure Q defined on J^t 
by 

:= exp{c(5T. ^ Ui) - -c^T\ ^ nj^} 



where 



c > max0-(O+) . 



The continuous process 

B[ :=Bt-ctJ2ni 
is a Q-Brownian motion on [0, T] and now 

dXt = dB[ - ^ nii/;'i{Xt.ni - ai)dt - ritdLt 

where 

ipi{u) := (f)i{u) — cu i eI . 

If Q{aj <T)=0 then F{aj < T) = and if this is true for any T we obtain F{aj = oo) = 1. 
c/ We are now going to prove that ai = ri = oo a.s. (with m > 2). Por any J C I let 



qj 



= span{nj,j G J} 

= dimVj (21) 
= orthogonal projection onto Vj . 



If qi = 1, then m = 2, ni + n2 = and Hi = Ki = 0. Assume now qi > 2 and Hi / 0. 
Choose some z e Hi and set 

Zt := TTiiXt - z) . (22) 
Then ^ ^ 

Zt = Zo + Ct - Y. [ ni(f>'i{Ui)ds + I HsdLs (23) 

where C is a gi-dimensional Brownian motion. Set 

St := iZtf . 

Then ^ ^ ^ 

St = So + 2[ Zs.dCs-2y2 [ UicP[{Ui)ds + 2 [ ZsMsdLs + qit. (24) 
Jo ^gj ./o ./o 

Prom Lemma 2 we deduce that on dD = Ujci^^j 

Zs.rzs = {Xs - z).ns = 
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and thus 

rt 



Let < T < oo. For t < rj A T, 

JO '-'s ^gj JO '^s Jo 

From the assumption made in b/ there exists < c < oo such that 0^ < on (0, c] and 

> l#:W)l<is (26) 

> — OO . 

We now proceed as in ([22], p. 47). As t ^ ri A T, the local martingale part in the r.h.s. of 
(25) either converges to a finite limit or oscillates between +oo and — oo. Thus it does not 
converge to — oo and a.s. Srj/\T > 0. Therefore 

P(ri < T) = 

and the conclusion follows since T is arbitrary. 

d/ Let now J C I with 2 < |J| < m — 1. We shall show by a backward induction on |J| 
that P(tj = oo) = 1. Remark that the backward induction assumption entails the equality 
a J = Tj a.s.. As previously done we may assume qj >2 and Kj ^ 0. Select now z £ Kj and 
set 

Zt := TTj{Xt-z) ^ ^ t ■ t (27) 

= Zo + Ct - J2jeJ lo (l^'j{Ui)ds - Y.i^3 /o ^J^* <P'AUl)ds + vrjn, dL, 

where C is a gj-dimensional Brownian motion. Let again St := \Zt\^. For e > and r > 
we set 

Ts := mf{t>0 : St+mim^jiUif < 2e^} 

Pr = inf{t >0 : \Xt\>r} . ^ ' 

From the induction assumption we infer that ^ oo as e goes to 0. Let < T < oo. We 
introduce the equivalent probability measure Q defined on J^t by 

_ 1 jr.^p^AT I ^^^^ 1^^.>^^<^^(C/]) ds} . 

Then 

Dt := Ct- V l{m>£}0-([7*)7rjnjds 

is a Qj-dimensional Q-Brownian motion on [0, T]. For t < A pr AT, 

St = So + 2 Zs.dDs - 2 /o Ul mi)ds - 2 Z,^, ^ui<e} Zs-n, mi) ds 
+2ElcI,L0J /o I^l(^s) Zs.UsdLs + qjt 

(30) 
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and for t < aj At^ A pr AT, 

log 5, = log5o + 2£^-2E,,./o*^^cl. 

+(^J-2)/of • 

From the induction hypothesis and the continuity of paths, if cij < oo for any L {Z! J there 
exists an interval (crj — S,aj] of positive length on which Xg Kj^. Therefore 

i-ajAT^AprAT V- „ 

- i^^(Xs)^^dLs>-(x>. (32) 

Jo "^s 

For s < Te, if {/] < e for some i J, then > and we obtain as well 

l{t/j<£} V Zs-Hids > -oo . (33) 
The other terms behave as in c/ and thus 

= Q{aj <T,AprAT)= P(C7J <TeAprAT) . (34) 

Letting e go to 0, r and T to oo we get 

¥{aj = oo) = 1 

and we are done. □ 



4 Keeping off from a wall 



We first recall some facts in the one-dimensional setting |21j . Let : M ^ (— oo,+oo] be 
a convex lower semicontinuous function. Assume 4> = +oo on (— oo,0) and on (0, +oo). 
Consider the one-dimensional equation 



dYt = dBt- (l)'{Yt)dt+\dL^t 
Yt > 

where is the local time of Y at 0. There are three types of boundary behavior: 



(35) 





repulsion 


(f){0) < oo 


weak: local time not zero 


(j){0) = oo, /Q^exp{2(/)} < oo 


middle: local time zero 


4>{0) = oo, jQ^exp{2cl)} = oo 


strong: boundary not hit 



We shall check the behavior of the multidimensional process X accords with this classi- 
fication in the neighborhood of the faces of the polyhedron. For any z G I we respectively 
write Hi,Ki,ai, n in place of i/{ j} , , cr{i} , r{ j} . 



Proposition 5 For any i G I such that (j)i{0) = oo and any t > 0, 

[ iH^XsJdLs = . 
Jo 



(36) 
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Proof. Prom the occupation times formula and Lemma 1 we obtain 

/•oo rt 

/ L-{Ui)Wi{a)\da= W,{Ul)\ds<^ (37) 

JQ Jo 

and from (j)i{0) = oo and the continuity of o h^- Lf{Ui) we deduce 

= . (38) 

Thus 

= ut-m+ 

= /q lH,{Xs)ni.dBs - /q lHi{Xs) Y,jei (f>'j(Ui) rii-rij ds + iRiiXg) rii-Us dLg 

= /q {Xs) Ui-Us dLs (39) 

= fo'^Ki{Xs)dLs 

= /o*l^,(X,)dL, . □ 
We now set for any i G I and x > 

/X 
exp{2(^i(u) - (f>i{l))}du . 

Theorem 6 For any i El such that Pi{0) = —oo or equivalently 

[ exp{20i} = oo, (40) 

then P(c7i = oo) = P(Tj = oo) = 1. 

Proof. Prom Ito formula and Proposition 5 we obtain 

PiiUi) = piiU^) + j'pf,{Ul)[dCl - rii.nj (k'M)'^' + E ^k,{Xs) n,.n^ dL,\ (41) 

where C* = B.Ui is a one-dimensional Brownian motion. As in the proof of Theorem 2, let 

T, := mi{t >0 ■.Ui + min,^,([//) < 2e} , . 

Pr = mi{t >0 : \Xt\>r} . ^ ^ 

Let < T < oo. We again introduce the equivalent probability measure Q defined on Tt by 

$ = exp{/;^^''"^^E,yil{t/|>,}<A;(^l)n..n,dC^ 

_ijre^Pr.T I ^^^^ l^^,^^^0;(t/i) n,.n,P d.} . 

Then 

Di:=Ci- Yl hui>e}^'jiUi) (44) 

is a Q-Brownian motion on [0, T] and for t < A A T, 

Pi{Ui)=pi{U^)+ /"^;(C/i)[di5^-^l|f;.<,}n,.n,>;.(C/i)ds + ^lK,(X.)n,.n,-dL,] . (45) 
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As in the proof of Theorem 2, for any j ^ i, 

r(Ti/\Ts/\pr/\T 

- ^{Ui<e} PiiUi) ni.rij (j)'j{Ui) ds > -oo (46) 

^0 

and 

+ iK.iXs) PiiUl) Ui.n, dL, > -oo (47) 

^0 

and then 

= Q(a, < A pr A T) = F{ai < r, A pr A T) (48) 
meaning that P((Tj = oo) = 1. □ 

5 Applications 

5.1 Brownian particles with nearest neighbor repulsion 

H.Rost and M.E.Vares [26j have considered the fohowing system: 

dXl = dBl + (j)' {Xj - Xl ) dt 

dXl = dBl + {(p'{Xl+^ - XI) - cP'iXl - Xf)) dt i = 2, . . . , n - 1 (49) 
dX't^ = dB^ - (p'iX]} - X^^^) dt 

where X^ < . . . < X" and cp is a positive convex function on (0, oo) satisfying 

</.(0) = oo, (/.(cx))=0, [\^'{x)fe-^^^''Ux<oo. (50) 

Jo 

This is a MSDE where function $ is given by (8) with (/)i{x) = <?!>(\/2 x), rii — — ^(gj+i — Si), 
ai = for i = 1, ... ,n — 1 and ej the j-th basis vector. Condition (50) for non-colhsion is 
stronger than (40) as can be seen from Schwarz inequahty: 

1 rl 

24, 



Jo Jo 



e 



5.2 Wishart and Laguerre processes 

Wishart processes have been introduced in [2] and [3j. If i? is a n x n Brownian matrix, 
a Wishart process with parameters n and 5 > n + 1 may be obtained as a solution to the 
matrix-valued SDE 

dSt = y^tdBt + dB't^/sl + 5Indt . (51) 
The eigenvalues process {X\, . . . , A") of {St] satisfies 



d\\ 



r- \* _(_ \i 

2^Jxi dWl + {5 + '^^ 1 < i < n , (52) 



and the square roots = y 



drl 
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where {W^, . . . , W") is a n-dimensional Brownian motion. N.Demni [T3] has remarked that 
this system is a MSDE with 

$(r\ . . . ,r"') = -^[(5 - n) ^log + ^log(r* + + ^log(r* - r-')] (54) 

i i>j i>j 

on {0 < < . . . < r"} and oo elsewhere. The system (53) has a strong solution for 6 > n. 
If 5 = n, we must add to the right hand side of (53) a local time at that disappears in (52). 
It has been proven in [3j that the eigenvalues never collide and if moreover 6 > n + 1 the 
smallest one never vanishes. This is in accordance with Theorem 6. 

Laguerre processes are Hermitian versions of Wishart processes. Only constants are 
changed in (52), (53) and (54). 



5.3 Reflection groups and Dunkl processes 

We only give a short introduction to this topic and refer to [19] and |25j for more details. For 
a G \ {0} we denote by Sa the orthogonal reflection with respect to the hyperplane Ha 
perpendicular to a: 

Ot.'X 

Saix) = X - 2—-TT . (55) 

A finite subset R C \ {0} is called a root system if for all a G R 

i? n Ra = {a,— a} ; , . 

Sa{R) = R. ^ ' 

The group W C 0{N) which is generated by the reflections {sq,, a € R} is called the reflection 
group associated with R. Each hyperplane Hp := {x G : p.x = 0} with /3 G \ Uagij-ffc 
separates the root system R into R+ and R-. Such a set i?+ is called a positive subsystem 
and defines the positive Weyl chamber C by 

C := {x G : a.x > Va G . (57) 

A subset S of R-^- is called simple if S" is a vector basis for span{R). The elements of S are 
called simple. Such a subset exists, is unique and we actually get 

C = {x G : a.x > Va G 5} . (58) 

A function k : R on the root system is called a multiplicity function if it is invariant 
under the natural action of W on R. If the multiplicity function k is positive on R^, we 
define the radial Dunkl process as the C-valued continuous path Markov process whose 
generator is given by 

C^uix) = lAn(x) + M«) (59) 
2 a.x 

for u G C^(C) with the boundary condition a.Vn(x) = for x G /^o. Then may be 
viewed as the solution to the MSDE 



dYt = dBt - V^{Yt)dt 
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where B is a Brownian motion and 



$(y) = ^(«) log(«-y) (60) 

on C and ^ = 00 elsewhere. It was proved in (|9j, ^lOj) that this equation has a unique strong 
solution and if moreover k{a) > 1/2 for any a € R then the process never hits the walls Ha 
of the Weyl chamber. In [14], it is proved that if k{a) < 1/2 for a simple root a, then the 
process hits Ha a.s. As a consequence of this result and of Theorem 6 (see also the statement 
at the bottom of p. 117 in [10]), we are in a position to classify the boundary behavior of the 
radial Dunkl process in the Weyl chamber. 

Proposition 7 For any a € i2+ let a a ■= inf{t > : G Ha}- 

• Ifa£R+\S, then F{aa = 00) = 1, 

• If a £ S and k{a) > 1/2, then P((Tq = 00) = 1, 

• If a £ S and k{a) < 1/2, then P((Tq, < cxd) = 1. 

5.4 Trigonometric and hyperbolic interactions 

Others interactions have been studied in [7]. 

The trigonometric system ([E], [18], [28]) reads 

dX{ = dB{ + ^ Em.- l<j<n (gi^ 

Xl <Xf < ...< <Xl + 2Tr 

This can be interpreted as the solution to the MSDE associated with 



where 



,/,(u) = -7 log ( sin ^) < -u < ^ ^g^^ 
= 00 elsewhere. 

It has been proved in [7] there exist a.s. collisions if 7 < 1/2. 
The hyperbolic system ([23], [27]) is 

dXi = dBl + 7 Efc^, coth (X,^' - X,fc) l<3<n 

xl<xi<...<xf. 



In this case 



with 



$(x) = E '^(^•^^) (65) 

l<jr<A:<n 



(^(n) = — 7 log (sinh(-v/2n)) u>0 

= 00 elsewhere. 

and collisions occur with positive probability if 7 < 1/2. 



(66) 
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